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A novel explanation of the relativistic self-induced transparency effect during superintense laser 
interaction with an overdense plasma is proposed. We studied it analytically and verified with 
direct modeling by both PIC and kinetic equation simulations. Based on this treatment, a method 
of ultrashort high-energy electron bunch generation with durations on a femtosecond time scale is 
also proposed and studied via numerical simulation. 
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Recent progress in laser technology has paved the way 
of exploring previously unattainable regimes of ultrain- 
tense laser-plasma interaction pQ. In up-to-date experi- 
ments, laser intensities of the order of 10 21 W/cm 2 and 
higher can now be achieved, implying that the goals like 
the fast ignition concept for inertial confinement fusion 
(ICF) PJ, laser-plasma based accelerators of charged par- 
ticles [3J 3] , and compact sources of short-pulse x-ray ra- 
diation [5] might soon be within reach. It seems that 
this is exactly where fundamental physics and technolog- 
ical progress intersect. Here we address the fundamental 
issue of ultraintense electromagnetic (EM) wave prop- 
agation through classically overdense plasmas, namely, 
why it occurs in the regime of relativistic self-induced 
transparency (SIT). In fact, this regime was first con- 
sidered in the pioneer work [B] in the form of stationary 
plane wave solution and later, in the 1970's, the station- 
ary solutions were extended to non-homogeneous plas- 
mas [7J 15] ■ Recently, another type of solution was found 
indicating that an ultraintense EM wave can penetrate 
into overdense plasmas over a finite length only, form- 
ing structured plasma distributions as a sequence of elec- 
tron layers separated by about half a wavelength wide 
depleted regions, so that this strongly nonlinear plasma 
structure acts as a distributed Bragg reflector 0. How- 
ever, such analytical studies, done in a rather academic 
manner, do not provide an answer of how and why the 
propagation of a laser pulse into overdense plasma oc- 
curs. Numerical simulations of ultraintense laser inter- 
action with overdense plasmas, on the one hand, have 
confirmed that the SIT effect takes place but also have re- 
vealed additional effects like anomalous longitudinal elec- 
tron heating, which can change the optical properties of 
the plasma and strongly influence pulse propagation dy- 
namics PH El DUU H US]- Nevertheless, the funda- 
mental question of the SIT regime - why the penetration 
into classically overdense plasma in the form of traveling 
wave occurs still has no reasonable explanation, espe- 
cially taking into account strong peaking in the electron 
density in the front of the laser due to the ponderomo- 



tive force pushing electrons forward [16 . It is generally 
recognized that the penetration occurs through lowering 
of the effective dielectric constant due to relativistic elec- 
tron mass correction and plasma heating up to relativistic 
temperatures under the action of laser field. 

In this Letter, in order to get an insight into the physics 
we focus special attention on the boundary electrons dy- 
namics and show that they are in unstable position, start- 
ing to be accelerated toward the incident wave. Thus, 
lowering of the effective plasma frequency occurs primar- 
ily due to the continuous generation of electron fluxes es- 
caping into vacuum and thereby decreasing the electron 
density, which allows the laser pulse to propagate fur- 
ther into the plasma. The subsequent dynamics of the 
escaped electrons strongly depends on their interaction 
with a slowly moving (with velocities ~ 0.1c, where c is 
the speed of light in vacuum) quasistanding wave struc- 
ture comprising the incident and reflected EM waves. 
Eventually they are divided into two groups of electrons. 
The first one is coming back into the plasma forming an 
electron beam with sufficiently high density. Duration 
of this electron beam can be very short if an ultrashort 
pump laser is used. The second group of electrons is 
captured into potential walls and, in fact, it can be rep- 
resented as two interpenetrating electron streams. Only 
these electrons gradually turn to heating. Indeed, the 
strong longitudinal heating takes place, especially for rel- 
atively high incident intensities, but it is rather concomi- 
tant effect than origin, although it can change even mode 
of the propagation. Based on the mechanism of electron 
beam generation we also propose here a new method of 
ultrashort electron bunch production on a femtosecond 
time scale. 

We shall start our analysis with the relativistic cold 
fluid model and consider a circularly polarized normally 
incident EM wave from vacuum (z > 0) onto a semi- 
infinite overdense plasma (z < 0). In a linearly polar- 
ized wave, inherent [j, B] heating at relativistic intensities 
(see, e.g., [17]) essentially complicates the classical effect 
of SIT or even makes it significantly different from the 
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FIG. 1: Schematic of quasistationary model of laser-plasma 
interaction: spacial distributions of electromagnetic field a, 
electron n and ion n; densities. 



latter. Since the relativistic SIT effect takes place on a 
time scale shorter than the ion response time, ion dynam- 
ics will be neglected. As was shown in [8j|9j[T8], there are 
two different cases with respect to the plasma density N a : 
(i) N < 1.5N C and (ii) N a > 1.5N C (N c = mu 2 /te 2 is 
the critical density and u> is the laser frequency), which 
are also confirmed by computer simulations |131 114] . In 
fact, at N < 1.5N C , as was shown by Marburger and 
Trooper [H], in a stationary regime there is a continu- 
ous family of solutions for the everywhere positive elec- 
tron density distribution, which means that penetration 
deep into plasma occurs in a "classical" way, through 
lowering the effective plasma frequency due to both rela- 
tivistic and ponderomotive nonlinearities. Modeling the 
SIT effect using PIC simulations and the fluid approach 
1 1 3] showed quantitatively close results, indicating that 
kinetic effects are not important for this case, at least 
for incident intensities not much exceeding the thresh- 
old of penetration [16]. The situation changes drastically 
at higher densities N > l.5N c , when the relativistic 
EM wave can considerably reduce the electron density 
in the region of its front until electron cavitation takes 
place. This electron cavitation acts as a "wall" making 
the plasma opaque to the propagation. In order to give 
an answer why the EM wave is able to propagate further 
deep into the plasma, we consider the dynamics of bound- 
ary electrons for the incident intensity near the threshold 
of the SIT regime. Our starting point is a solution given 
by Cattani et at |16j . which can be depicted schemat- 
ically as in Fig. [T] Let us obtain the electron density 
distribution from the Poisson equation using the balance 
of ponderomotive and electrostatic forces ((f)' = 7'): 
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Here we introduced normalized variables: 
Re[a(e x + ie y )e luJt ] is the vector potential; <fr = eip/mc 2 
the scalar potential; n — N e /N a the electron density; 
n o = N a /N c the overdense parameter; 7 = \f\ + a 2 the 
relativistic factor, the primes denote differentiation with 
respect to the longitudinal coordinate £ = luz/c. By us- 



ing the equation for vector potential, as in Ref. [IB] , 
for the incident intensity equal to the threshold value 
we arrive at the following equation for the field at the 
boundary where electrons vanish (a(£b = £t>*) = a&*): 



2(n c + a b J - 3n Q Wl + 



0. 



(2) 



Let us now show that the electrons at this boundary 
are in unstable equilibrium state. To do this, we take the 
ponderomotive force acting on them in the form 



F v = -i = - 
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(3) 



and calculate its gain when the electrons arc shifted by 
a small distance A£ opposite to the propagation direc- 
tion, i.e., to the vacuum region (£ > Then for EM 
field we can use the vacuum solution, which represents 
a sum of incident and reflected waves, that should be 
self-consistently matched at the boundary with the field 
inside the plasma. Substituting this solution into Eq. ^ 
and taking small variation from both parts, we obtain 



AF„ = 2 



(4) 



where aj, is the value of the field at the electron layer 
boundary. For the electron to be in stable state when it 
is shifted from the equilibrium position, the total force 
acting on it must return it back to the equilibrium po- 
sition. For that, the ponderomotive force gain must be 
more than the electrostatic force gain: 



AF p > -n A£. 



(5) 



From this we find the following condition on the field at 
the boundary: 



l + a 2 b > 2(n + al). 



(6) 



Comparison of this expression with Eq. Q shows that 
the stability condition ^ ceases to be valid for the value 
of the field at a& = a^, i.e., for the incident intensity 
equal to the threshold. This conclusion is highly impor- 
tant as it shows that all the skin-layer solutions are not 
only equilibrium (as follows from the balance of forces), 
but stable too. When the intensity of incident radia- 
tion exceeds the threshold, the electrons at the boundary 
become unstable and start to move toward the incident 
wave. This motion does not cease at least until electrons 
escape from the plasma region (0 > £ > £&), which can 
be seen, for example, from Eq. Q. This occurs because 
the gain in ponderomotive force only decreases with in- 
creasing field and, in spite of the fact that there may 
exist intervals where the field is less than at the electron 
layer boundary (near the standing wave nodes), this is 
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insufficient for the balance of force. The balance is pos- 
sible only when electrostatic force changes its sign, i.e., 
outside the plasma layer, £ > 0. 

Direct calculation of ponderomotive force gain at the 
electron layer boundary with a small increase of incident 
intensity may be an additional argument confirming that 
the electrons start to move towards the laser, when the 
threshold is exceeded. Indeed, the substitution a — > a + 
ai into Eq. where |ai| « \a\, yields the following 
gain at the threshold intensity: 



2n (y/l + al-l)-al 

As seen from Eq. ([7]), the gain is negative, i.e., the pon- 
deromotive force decreases and becomes less than the 
electrostatic force. Thus, with the excess over the thresh- 
old, the electrons start to escape from the plasma to vac- 
uum and may get even far outside. The maximal energy 
that may be acquired by the electrons is proportional 
to the magnitude of the charge separation field, which 
grows with increasing n , in particular, at n 3> 1 it can 
be estimated as E m 1.3nl [TB"] . 

To consider full scenarios of the SIT propagation in- 
cluding kinetic effects, we have done a number of fully 
relativistic ID PIC simulations for a wide range of plasma 
densities, n = 1 — 100, paying special attention to the 
boundary electron dynamics. The incident laser pulse 
was taken either Gaussian or in the form of a unit 
function with smoothly growing amplitude up to max- 
imal value exceeding the threshold. First, simulations 
showed good agreement with the SIT threshold values 
given in |16j . which clearly indicates the validity of the 
cold fluid plasma model at intensities below the thresh- 
old. For these intensities the electron density at the 
vacuum-plasma boundary is redistributed under the ac- 
tion of ponderomotive force and a skin-layer structure 
is formed, as in Fig. [T] They are similar to the field 
and electron density distributions of the dynamic prob- 
lem presented in Figs. |2ja), [3ja), until the field at the 
boundary exceeds the threshold. In Figs. 2 and 3 the 
dynamics of laser propagation for two values of plasma 
density n — 2 and n a — 5 is shown, respectively. Here 
analysis is given for a semi-bounded incident pulse with 
rising time equal to 3 laser periods (time and space vari- 
able are normalized as: t — > uit, z — > uz/c) and with 
maximal amplitude (ao = 1.7 and 12.8) slightly exceed- 
ing the threshold (a t h = 1-58 and 12.6, respectively). At 
the beginning, simulations show that the dynamics of the 
interaction is quasi-stationary, i.e. the electromagnetic 
field structure having the form of a quasi-standing wave 
with a skin-layer distribution in plasma gradually pushes 
electrons deep inside the medium, creating an uncom- 
pensated ion charge at the plasma boundary and thus 
increasing the longitudinal quasistatic field. When the 
incident intensity becomes higher than the threshold, vi- 



olation of the quasi-stationary interaction starts, in the 
first place, from the acceleration of boundary electrons in 
the skin layer towards the incident wave, which is read- 
ily traced in the phase portrait of longitudinal motions 
(see Figs.|2jb,d,f,hj) and[3]jb,d,f,h)) and is in qualitative 
agreement with the above analytical considerations. The 
difference between Fig. [2] and Fig. [3] is that for higher 
plasma density the electrons acquire higher energies and 
may propagate at large distances from plasma into vac- 
uum. As the electrons escape, the plasma eventually ac- 
quires a positive charge, thus the electrons get into the 
total field that is a sum of the quasi-static field and the 
standing wave which can be regarded as a sequence of 
potential wells formed by the ponderomotive potential. 
With a small excess over the threshold, the electrons 
both oscillate, primarily in the first potential well, and 
return to plasma forming electron beams moving deep 
into plasma as is shown in Figs, g^hj) and gf.h). It 
should be noted that subsequent dynamics of the escaped 
electrons is quite complex and strongly depends on laser 
and plasma parameters. However, eventually there are 
two groups of electrons. The first one is coming back 
into the plasma in the form of an electron beam that can 
propagate deep into the plasma and has sufficiently high 
density. Moreover, the duration of this electron beam 
may be very short if an ultrashort pump laser pulse is 
used. The second group of electrons is captured into po- 
tential walls and, in fact, it may be represented as two 
interpenetrative electron streams, clearly indicating that 
the hydrodynamic approach is violated. Only these elec- 
trons gradually turn to heating. This is clearly seen on 
the right-hand side of Figs. [2] and [3j where the phase tra- 
jectories inside the separatrices corresponding to poten- 
tial wells are gradually densely filled everywhere, which 
is the evidence of effective electron heating up to the 
temperature corresponding to the ponderomotive poten- 
tial. With a substantial excess over the threshold, the 
electrons may move over large areas of a standing wave. 
Note that the presence of oscillatory trajectories for con- 
tinuously incoming electrons in ponderomotive potential 
is indicative of the need to use the kinetic description be- 
cause the electron fluxes moving in forward and backward 
directions intersect here. It is also worth noting that the 
electron beam moving deep into plasma has a quite high 
density, only an order of magnitude less than the initial 
plasma density and energy of several MeV; parameters 
of the beams corresponding to Figs. [2] and [3] are the fol- 
lowing: density m ~ 0.2 and 1.3, maximal kinetic energy 
(in units of the rest energy mc 2 ) tbmax — 0.63 and 19.3 
with energy spread of Ae;, = 0.4 and 1.3, respectively. 

The process of acceleration of electrons from plasma 
to vacuum and then back is evidently accompanied by si- 
multaneous wave penetration deeper into plasma. In this 
case, the wave reflects from a plane gradually moving into 
a plasma that separates the transparent and nontrans- 
parent regions. The penetration velocity of the incident 
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FIG. 2: Space-time dynamics (a,c,e,g,i) of the interaction of 
incident (on the right of the boundary) semi-bounded pulse 
of maximal amplitude ao = 1.7 with supercritical plasma 
no — 2: distribution of electron density (continuous curve), 
ions (dotted curve), absolute value of vector potential (dashed 
curve), scalar potential (dot-dash curve) at different moments 
of time. The corresponding phase portraits of the electrons 
are plotted on the right (b,d,f,h,j) 



FIG. 4: The same as in Fig. [2] but for the incident 5-fs Gaus- 
sian pulse with ao = 60. Data for the generated electron 
beam: (g) snapshot of the density distribution and (h) its 
detailed phaze portrait 
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FIG. 3: The same as in Fig. [2] but for ao = 12.8 and n Q = 5 



wave, or to be more exact of the reflection plane, depends 
on the excess over the threshold of SIT. If the threshold 
is exceeded substantially, the velocity is v « 0.1 — 0.2c, 
which agrees with data in jTUl H21 US]- However, for a 
small excess over the threshold, after penetration to a fi- 
nite depth, the interaction enters a quasi-stationary stage 
close to that described in the hydrodynamic approxi- 
mation |13j , except for the absence of purely cavitation 
regions. These regions are partially filled by electrons, 
which indicates at least that for high local ponderomo- 
tive potential the effective electron temperature is high. 
In view of the above result of electron beam accel- 



eration, we propose here a new method of high-energy 
electron beam generation with ultrashort durations. The 
first advantage of this method is that the density of the 
beam can be very high, for example, for the parameters 
as in Fig. [2] it is of the order of 0.1n o , just ten times less 
than the solid density but with energy of several MeV. 
On the other hand, by using superintense ultrashort laser 
pulses, electron bunches with a few femtosecond or even 
shorter duration may be produced. The interaction of a 
5-fs laser pulse of intensity I ~ 10 22 Wcm~ 2 (a = 60) 
with plasma of n Q — 2 is depicted in Fig. [4j where an 
electron bunch generated at the plasma boundary has a 
sub-fs duration, energy of order GeV and maximal beam 
density of 10 22 cm~ 3 (rib ~ 10). Results of simulations 
for various incident intensities are systematized in Ta- 
ble 1 where parameters of the generated electron beams 
are listed: nb max - maximal electron density, < rib > - 
beam density averaged over the time interval indicated 
in the last column, the same for the kinetic energy of 
the beam tb- Two averaging times are used because the 
temporal structure of the generated beams exhibits quite 
a narrow peak and a relatively extended pedestal, as is 
shown in Fig.jijg); it allows us to present results on elec- 
tron bunch structure in a better way. 
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In conclusion, we have presented a new treatment of 
the relativistic self-induced transparency effect that al- 
lows ultraintense laser pulses to penetrate deep into over- 
dense plasmas. This treatment is based on the careful 
consideration of the boundary electrons dynamics, which 
for incident intensities exceeding the threshold start to 
move towards the laser, thus decreasing the electron den- 
sity in front of the pulse. The subsequent dynamics of the 
generated electron beams strongly depends on the laser 
and plasma parameters, in particular, they may be con- 
sidered as ultrashort high-energy electron bunches. In 
view of such electron bunch dynamics, we also propose 
here a method of high-energy beam generation in a fem- 
tosecond time scale durations or even shorter. 
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